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Abstract 

We consider primordial fluctuations in thermal inflation scenario. Since the ther- 
mal inflation drives about 10 e-folds after the standard inflation, the time of horizon- 
exit during inflation corresponding to the present observational scale shifts toward 
the end of inflation. It generally makes the primordial power spectrum more deviated 
from a scale-invariant one and hence renders some models inconsistent with observa- 
tions. We present a mechanism of generating the primordial curvature perturbation 
at the end of thermal inflation utilizing a fluctuating coupling of a flaton field with 
the fields in thermal bath. We show that, by adopting the mechanism, some infla- 
tion models can be liberated even in the presence of the thermal inflation. We also 
discuss non-Gaussianity in the mechanism and show that large non-Gaussianity can 
be generated in this scenario. 



1 Introduction 



The thermal inflation [1-7], a mini-inflation which occurs long after the standard infla- 
tion, has been discussed, in particular, as a possible solution to the cosmological moduli 
problem [8-10], in which moduli particles would lead to various cosmological difficulties 
such as destroying light element synthesized by BBN, too much contribution to X(7)-ray 
background radiation, overclosure of the universe and so on. The thermal inflation can 
produce very large entropy to dilute the moduli density sufficiently and thus can be a 
solution for the moduli problem. 

Here we would like to discuss another aspect of the thermal inflation. When one con- 
siders the thermal inflation, primordial fluctuations are usually assumed to be generated 
from the primordial inflation. However notice that the thermal inflation drives the number 
of e-folds of about 10 due to the large entropy production. This means that fluctuations 
corresponding to the present observational scales exit the horizon at later time compared 
with the case of no thermal inflation. Since the time of horizon-exit becomes closer to the 
end of inflation, the primordial power spectrum is rendered to be more tilted in most cases. 
Hence even an inflation model consistent with observations in the absence of the thermal 
inflation can be pushed to outside the allowed range. For example, here let us consider 
the quadratic chaotic inflation model, where the inflaton potential is given hj V (x (p'^. 
In the usual case, the curvature perturbations probed by cosmological observations are 
assumed to correspond to the fluctuations which exit the horizon when the number of 
e-folds is A^inf = 50 — 60. By adopting this number, the spectral index is rig ~ 0.96 and 
the tensor-to-scalar ratio is r ~ 0.13, which is consistent with current observations such as 
WMAP5 [11]. However, when the thermal inflation occurs after the standard primordial 
inflation, the number of e-folds reduces by about 10, and then the spectral index is more 
red-tilted and the tensor-to-scalar ratio is more increased. When A^inf = 40, the tensor- 
to-scalar ratio is r ~ 0.2, which is inconsistent with current observations and invalidates 
the model. Thus in this respect, the thermal inflation can also affect the predictions of 
primordial fluctuations. 

In fact, the primordial fluctuations are not necessarily generated from the inflaton 
fluctuations. Another light scalar field such as the curvaton [12-14], modulus in the 
modulated reheating scenario [15,16] and so on are also known to generate (almost) scale- 
invariant and adiabatic primordial fluctuations consistent with observation^*^- Although 
by adopting these scenarios, one can alleviate the above mentioned issue, however, here 
we propose another mechanism which can naturally arise in the framework of the thermal 
inflation. During the thermal inflation, the effective potential of a flaton, a flat direction 
in super symmetric models, is lifted up by a thermal effect due to the couphng between 
the flaton and particles in thermal bath. The end of the thermal inflation is controlled by 
the strength of the thermal effect, or the coupling. If the coupling depends on some other 
scalar field and this scalar field fluctuates, the end of the thermal inflation also differ from 
place to place in the Universe via the fluctuations of the coupling. If the scalar field is 

Other mechanisms have also been discussed in the hteratures [17-20]. 
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light enough during inflation, almost scale-invariant density fluctuations can be achieved 



similarly to the case of the modulated reheating^ 



Another important issue in the primordial fluctuations is the Gaussian nature of the 
fluctuations. In fact, although recent observations are consistent with purely Gaussian 
fluctuations, but they may suggest that primordial local type non-Gaussianity^i is large 
compared to that of the standard (single-fleld) inflation model [11,21]. Since some of the 
alternative mechanisms such as the curvaton and the modulated reheating scenarios can 
also produce large non-Gaussianity, they have been the subject of intense study [22,23]. 
We show that large non-Gaussianity can be generated in the mechanism as well. Thus the 
mechanism proposed here would be interesting in this respect too. 

The structure of this paper is as follows. In the next section, we briefly describe 
the thermal inflation model which we consider here. Then in Section [3l we summarize 
the formalism to discuss the curvature perturbation in the model and some observational 
quantities such as power spectrum, bispectrum and trispectrum or nonlinearity parame- 
ters. In Section HJ we present a mechanism of generating the curvature perturbation at 
the end of thermal inflation and also discuss its non-Gaussianity. In Section [5], we discuss 
the current observational constraints for several simple inflation models in our scenario. 
The flnal section is devoted to summary of this paper. 



2 Basic Picture of thermal inflation 

In this section, we briefly review the idea of the thermal inflation [5] . The thermal inflation 
can be realized by utilizing a flat direction which exists in supersymmetric theory. Let us 
call such a direction as a flaton fleld and assume that the potential of the flaton fleld is 
given by 

= vo - i™v + ^ . (1) 

where the cfP' term comes from soft supersymmetry breaking and we have neglected the 
higher order nonrenormalizable terms (0^""*"^, n > 2). The VEV of is given by 0vev = 
■i^fiaton = A~"'^/'^m^/^Mp(^ and Vq = ""^^"yflaton/S to have V^(f flaton) = 0. When the flaton 
fleld interacts with some particles in thermal bath, such an interaction gives a thermal 
contribution to the potential of the flaton fleld as 

Vt = |t^0' , (2) 

In Ref. [24] the author has presented a mechanism of generating large scale curvature perturbation 
through the inhomogeneous cosmological phase transitions in the early universe. The mechanism proposed 
here is similar to this kind. 

Two types of non-Gaussianity are often discussed in the literatures, one is the so-called local type 
and the other is the equilateral type. In this paper, we only focus on the local type non-Gaussianity. 
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where g is an effective coupling between the ffaton and particles in thermal bath and T is 
the cosmic temperature. Then the total effective potential of the flaton field in thermal 
bath is given by 

y^^ = Vo + l{9T'-m^)<t>' + ^^f. (3) 

After the standard primordial inflation ends followed by the reheating, the cosmic temper- 
ature decreases as the Universe expands. Then, at some temperature when Vq gets larger 
than the background radiation energy density, the Universe is dominated by the false vac- 
uum of the flaton's potential, which drives a mini-inflation. When the temperature of the 
Universe has dropped down toT = Tc = m/ ^/g, the flaton starts rolling down to the VEV 
and then the mini-inflation ends. The duration of this thermal inflation can be expressed 
by using the e-folding number as 

iV(t„ ti„) = Hdt = In (^^^ = - In (^^^ , (4) 

where H denotes a Hubble parameter. Here Tin is the temperature at the onset of the ther- 
mal inflation, which is defined as Pr(^in) = Vo = "^^^Haton/^! i-e., Tin ~ "^^''^^filton' where 
Pr represents the radiation energy density. Taking some typical values of m ~ 10^ GeV 
and fflaton — 10^° GeV, we have Tin — 10^ GeV and Tc ~ 10^ GeV and hence the size 
of the e-folds during thermal inflation is estimated as Nit^taT) ~ 10, which corresponds 
to the number AA'th to be subtracted from the standard one. More precisely AA'^th is 
calculated by considering entropy production due to the thermal inflation. The total en- 
tropy S of the Universe after thermal inflation increases by a factor As = s{Tr)af/s(Tc)al 
where and Oc are the scale factor at the reheating after the thermal inflation and the 
end of the thermal inflation, respectively, s is the entropy density and T^ is the reheating 
temperature after thermal inflation. Since the flaton behaves like a matter after the end 
of the thermal inflation, the energy density of the flaton just before the reheating is given 
by Pflaton(Tr) = Vo{ac/ar)^. This gives the radiation energy density just after the reheat- 
ing, hence the entropy density can be written as s{Tr) = 4:pr{Tr)/3Tr. Then, A5 can be 
estimated as 

where 100) counts the effective degrees of radiation. Then AA^th is estimated as 

3 6N formalism and observational quantities 

Based on 6N formalism [25,26], the curvature perturbation on the uniform energy density 
hypersurface, (, on super-horizon scales is given by the difference of the neighbor back- 
ground trajectories, which is parametrized by the e-folding number, measured between 
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the initial flat hypersurface and the final uniform energy density hypersurface. That is, 
we have 

C(t/)^5iV(t;,t,) , (7) 

where N{tf, t*) denotes the e-folding number measured between those aXt = tf and t = t*. 

When the curvature perturbation ( originates from fiuctuations of a single scalar field 
a, ( at t = tf is given, up to the third order, by 

C{tf) = NM* + \N,,{6a,f + ^N,,,{6a,)\ (8) 

where N^j = dN/da^, and so on. The power spectrum Pq, the bispectrum B(^, and the 
trispectrum are given by 

{C^,^C^,^) = i27rfP^{h)6{h + h), (9) 

iCkSk^CkJ = W'Bdh, h)5{h + k2 + h), (10) 

iCkSk^Ck^Ck,) = i2nj%ih, k2, h, h)S{h + h + h + h), (11) 
where 5^ and can be expressed with the power spectrum as 

Bcih, k2, h) = ^/nl (Pc(^i)n(^2) + Pc(^2)Pc(^3) + Pc{h)P^{h)) , (12) 

T^(/ci, /C2, ^3, ki) = tnl (Pc(^13)-Pc(^3)-Pc(^4) + 11 perms.) 
54 

+^^7NL {Pc{k2)Pdh)Pdh) + 3 perms.) . (13) 

Here /nl, tnl and r7NL are the non-linearity parameters often discussed in the literatures. 
By adopting SN formalism, these parameters can be calculated as [27] 

^NL-^, ^^^NL-^. (14) 

When the curvature perturbation is generated from a single source, the relation tnl = 
(36/25) /^L holds, which is the case in the following discussion. Thus when we investigate 
the non-hnearity of fiuctuations, we do not consider tnl, but discuss /nl and (7nl- 

4 Curvature perturbation generated at the end of 
thermal inflation 

Now we discuss the generation of density fiuctuations at the end of thermal infiation. We 
assume that the couphng g depends on a light scalar field a, namely, g = g{cr). Then 
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the coupling g can fluctuate due to the fluctuation of the light field a which originates to 
quantum fluctuations during inflation and hence the fluctuation of the coupling g gives rise 
to the inhomogeneous end of thermal inflation. In such a case, from Eq. (jlj) the curvature 
perturbation from fluctuations of a at the end of thermal inflation can be given by, up to 
the third order. 



C=5N: 




5al 
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5al 



(15) 



where a prime denotes the derivative with respect to a* 
Now the power spectrum is written as 



Here we have used Tc oc g 



-1/2 



2^ 



(16) 



Since a does not evolve much during inflation, the scale dependence of the power spectrum 
comes from the time variation of the Hubble parameter during inflation. Then the spectral 
index is given by 
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-2e 



(17) 



where e is a slow-roll parameter and defined by e = (Mp[/2)(V0/y)^. Here V is the 
potential for the inflaton </> and = dV/dcf)^. Notice that in the standard inflation 
scenario, the spectral index is given by ra^ — 1 = — 6e + 2rj with 77 = M'^^iV^^/V) being 
another slow-roll parameter. 

To discuss non-Gaussianity of the curvature perturbation in the mechanism, we give 
the non-linearity parameters defined in the previous section. For the bispectrum, one 
usually uses the parameter /nl which is given by 



6. 

7/nL 
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9"l9 



V{9'/9r 

For the trispectrum, the non-linearity parameter (/nl is 



54 
25 



fi'NL 



9" /9 o 9/9 



(19) 



To investigate observational quantities in more detail, we assume a simple functional 
form of g as 



9 = 9o\^ 



2M2 



(20) 



5 



with go and M being a coupling constant and some scale, respectively. With this form of 
g, the curvature perturbation generated from the fluctuations of a is given by 



In general, however, fluctuations from the inflaton can also be generated. The curvature 
perturbation d^i from the inflaton can be expressed as 

Thus the ratio of the amplitude of the power spectra between those from the mechanism 
proposed here and the inflaton is given by 



PCi„f 2 \MJ \ M 



(23) 



On the other hand, if we assume that {a^^/M)^ <^ 1, the non-linearity parameter /nl can 
be written as 

The constraint on the non-linearity parameter from the WMAP5 [11,21] is /nl < O(IOO), 
which indicates that a./M > 0(10"^). By using above equations, we can rewrite the ratio 
in Eq. ([23]) as 



P 



Hence as far as e is not extremely small, the fluctuations from a in this mechanism can 
be responsible for cosmic density perturbations today. 

When the fluctuations from a are fully responsible for the scalar primordial fluctua- 
tions, the tensor-to-scalar ratio is 



MJ * VlOO/ V IQiOGeV 



(26) 



Thus, as other mechanisms of generating primordial fluctuations such as the curvaton 
and modulated reheating models, the tensor-to-scalar ratio is generally very small in this 
mechanism. 

Here we make some comments on (^nl- With the form of g assumed in Eq. (l20l) and 
when g'" is negligible, we have the relation between /nl and g-^i,: 

10 50 

5'NL — — ^/nl - ^ • (27) 



Thus in this model, the sizes of /nl and (^nl are almost the same. 

In the above discussion, we have neglected the evolution of the light field a after 
the horizon crossing until the end of thermal inflation. This assumption is valid when 
rrtfj <C i^th with iJth being the Hubble parameter during the thermal inflation. As we 
discussed in Section El ifth can be estimated as ifth = V^^"^ /M-p\ ~ 10^ eV and hence a 
has to be light as nifj ^ 10'^ eV. Moreover, a gets the effective mass from the coupling 
~ g{a)(fP'T'^. Thus such an effective mass should be also smaller than the Hubble parameter 
during the thermal inflation. In our scenario, however, this effective mass coming from 
the coupling can be negligible until the end of thermal inflation because the flaton field is 
trapped at = during thermal inflation phase. 

5 Primordial fluctuations in scenarios with thermal 
inflation and models of inflation 

Now in this section, we discuss what inflation models can be compatible with the thermal 
inflation in the light of current cosmological observations. As mentioned in the intro- 
duction, when a mini-inflation occurs after the primordial inflation, the time when the 
fluctuations of the reference scale in the present observations exit the horizon during in- 
flation is shifted toward the end of inflation. This, in most cases, means that the spectral 
index is more tilted. Thus even if an inflation model consistent with current observations 
in the standard case can be excluded when the thermal inflation occurs. However, this 
situation can be relaxed by adopting the inhomogeneous end of thermal inflation discussed 
in Section H] since fluctuations from the mechanism are generally more scale-invariant as 
seen from Eq. (ITTj) . We investigate this issue by working out explicitly some inflation 
models and discuss in what cases they are consistent with current observations. 

5.1 The number of e-folds 

To predict the quantities such as the spectral index and tensor-to-scalar ratio, we need to 
specify the number of e-folds when the reference scale exits the horizon, at which above 
mentioned quantities are measured. The number of e-folds at the time when fluctuations 
of the scale k exit the horizon can be given by 



the energy density of radiation at the reheating after the primordial inflation. Here we 
have assumed that the inflaton oscillates around a minimum of a quadratic potential for 
some time after the inflation ends. 



where V- 
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When the thermal inflation occurs, A'^inf given above would be reduced by AA'^th of 
Eq. (Q. Thus, in this case, the number of the e-folds is modified as 



/ A: \ 2, f W^GeV\ 1, Ao^GeVX 

Since the spectral index and the tensor-to-scalar ratio depend on the number of the e- 
folds, some models with the reduced Nij^f due to the thermal inflation may be excluded by 
current cosmological observations. 

In the following, we consider some representative models of inflation such as chaotic 
inflation, new inflation and hybrid inflation models. First we focus on how the spectral 
index and tensor-to-scalar ratio are modified in the presence of the thermal inflation when 
primordial fluctuations are generated from the inflaton and discuss the compatibility with 
observational constraints. Then we investigate the case where primordial fluctuations are 
produced by the inhomogeneous end of thermal inflation, which is discussed in Section IU 
and see that some excluded inflation models can be relaxed by using this mechanism. 



5.2 Chaotic inflation 

Chaotic inflation models have the potential of the following form: 



VW = AMj^^j . (30) 

where (j) is the inflaton field and n is assumed to be even integer to have a minimum at 
the origin. The spectral index and the tensor-to-scalar ratio in chaotic inflation models 
are respectively given by, in terms of iVinf , 

n, — 1 = , r = . 31 

From the above expressions we can find that as A^inf becomes smaller, Ug is more red-tilted 
and r becomes larger. In Fig. [H contours are shown for 68 % (green dashed line) and 95 % 
(red solid line) CL. allowed region derived from WMAP-I-BAO-I-SN. The blue dotted line 
corresponds to the predictions of chaotic inflation models on the n^-r plane for n = 2 (left) 
and n = 4 (right) cases, respectively. The values of Ug and r in the standard inflation 
case are also shown with small circles for A^inf = 30,40,50 and 60. When A^i^f < 45, we 
find that the quadratic model is excluded at 95 % CL. For the case of n = 4, as is well- 
known, the model is ruled out even if A^inf = 60. However, if fluctuations come from the 
inhomogeneous end of thermal inflation, the situation drastically alters. As discussed in 
the previous section, the spectral index is given by — 1 = — 2e and the tensor-to-scalar 



8 



ratio would be negligibly small. In Fig. [H we also plot the predictions of and r from 
the inhomogeneous end of thermal inflation with small squares. Although the quadratic 
chaotic inflation model with iVinf < 45 are excluded by current observations when the 
infiaton generates primordial fluctuations, such a model becomes viable in the present 
mechanism even if A^inf is small as 30. More interestingly, the quartic inflation model, 
which is excluded with its orig inal form, can be viable in the present mechanism as shown 



in the right panel of Fig. [j^ 
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Figure 1: Contours indicate 68 % (green dashed line) and 95 % (red solid line) CL. 
allowed regions derived from WMAP+BAO+SN [11]. The blue dotted line corresponds to 
the predictions for the quadratic (left) and the quartic (right) chaotic models. The small 
purple circles and red squares represent the values of Ug and r from the standard inflation 
and inhomogeneous end of thermal inflation, respectively, for Ajnf = 30, 40, 50 and 60. 



5.3 new inflation 

For the new inflation model, we consider the following potential: 



V = Va 



l-2( ^ 

V 



2n 



(32) 



where v denotes the VEV of the infiaton. We assume that (l)/v ^ 1 during inflation. 
Then, the slow-roll parameters are given by 



e = 2n' 



2(n-l) 2 



7] = — 2n(n — 1) 



(33) 



For a similar argument in the framework of the mixed infiaton and curvaton or modulated reheating 
scenario, see [28]. 



9 



From these expressions, we find that ^/rj = <^ 1. Hence the tensor-to-scalar ratio 

is very small and the spectral index is given by 

ra, - 1 ~ 2r/ . (34) 

The e-folding number and the value of the inflaton field at the time of horizon crossing 
are related by 



r 1 



V 



inf 



4 VMpi 
1 

, 2n{n — 



In ^ 



1 



(n = 2) 



M, 



pi 



n 



n 



{n > 3) 



From these equations, the spectral index is given by, in terms of A'^inf, 

[n = 2) 

[n > 3) 



inf 



n 



n 



2 N, 



inf 



From the above expression, we find that when the curvature perturbation is generated 
from the inflaton, Irig — 1| becomes larger as A'lnf becomes smaller. For examples, when 
77, = 3, 77-5 = 0.933, 0.92 and 0.9 for the cases with A'^inf = 60, 50 and 40, respectively. 
From WMAP Syr data, the hmit on Us is obtained as 0.925 < ra^ < 0.997 (95 % C.L.) for 
negligible tensor-to-scalar ratio [11]. Hence, the models with A^i^f = 50 and 40 lie out of 
the 95 % C.L. allowed region. 

Now we consider the primordial fluctuations from the inhomogeneous end of thermal 
inflation. As discussed in the previous section, the condition -Pc„/-P<^i„f ^ 1 should be 
satisfied when the fluctuations from the mechanism dominate over those from the inflaton. 
However, this requires e(o"^,/M)^(Mpi/M)^ ^ 1, thus when e is very small as in the case of 
the new inflation, the above condition is difficult to be satisfied. Even if the fluctuations 
from the inflaton give a negligible contribution and those from the inhomogeneous end 
of thermal inflation dominate, its power spectrum becomes fairly scale-invariant since the 
spectral index is given by — 1 = — 2e <^ 0(0.01). Since the tensor-to-scalar ratio is also 
very small, such fluctuations are also excluded. Thus, the new inflation may not be well 
suited to the framework of the thermal inflation. 



5.4 Hybrid inflation 

Next we consider hybrid inflation. Although there are some variations for the potential 
for the hybrid inflation, we discuss some representative ones here. First we consider the 
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potential of the following form: 



V{4>) = K) + (35) 



When Vo ^ ^\(t>^i the slow-roll parameter can be written as 

e~— r/~^-^. (36) 



Since the condition Vq ^ ^^1^0^ implies that e <^ 77, the spectral index of this model can 
be written as — 1 = 2?7, which is blue-tilted since rj is positive here. The smallness of e 
also indicates that the tensor-to-scalar ratio is also negligible in this model, and in such a 
case the blue-tilted spectrum conflicts with the observations (see Fig. [1]). Thus, this type 
of hybrid inflation is excluded by current observations in its original form. As in the case 
of the new inflation, since e is very small in this model as well, it could be difficult to have 
the case where fluctuations from the inhomogeneous end of thermal inflation dominate 
over those from the inflaton. Furthermore even if the fluctuations from the mechanism 
dominates over those from the inflaton, the spectral index is 2^ 1, thus this again 
contradicts with current observations. 

Another possible realization of the hybrid inflation is so-called mutated hybrid inflation 
[29], in which the potential can be written as 

V{<P) = V,-^,Ml,(^^y , (37) 

where /i is a constant. The potential of this kind can appear, for example, smooth hybrid 
inflation [30] for n = 4. With this potential, the slow-roll parameters are 

where we assumed that Vq ^ I^M^^ ( • '^'^^^ condition also implies that e <S |?7|, thus 

the spectral index can be given by 

n.-l^-Z^-y) . (39) 

When n = 2, the spectral index is = 0.95,0.962,0.97 and 0.975 for A^'i^f = 30,40,50 
and 60, respectively. Thus even if the number of e-folds is reduced due to the thermal 
inflation, this model is not excluded. In fact, the central values of the allowed range for 
the spectral index from WMAP5 analysis is Us ~ 0.96, thus a slight reduction of A^inf may 
be favored. It should be noticed that the slow-roll parameter e in this model could also 
be very small, the fluctuations from the inhomogeneous end of thermal inflation may not 



11 



give a dominant contribution for the curvature fluctuations. Even if fluctuations from the 
inhomogeneous end of thermal inflation are the main component of (, the spectral index 
is again too scale-invariant to be consistent with observations, which was also the case for 
the new inflation and hybrid inflation with the potential of Eq. (l35l) . The mutated hybrid 
inflation discussed here can be compatible to cosmological observations in the framework 
of thermal inflation even in its original form. 

6 Summary 

In this paper, we have presented a mechanism of generating primordial curvature pertur- 
bation at the end of thermal inflation by considering the case where the coupling of a 
flaton field with the fields in thermal bath can fluctuate. We also show that there is a 
possibility of generating large non-Gaussianity in this scenario. 

We have also investigated the constraint on inflation models in the case where the 
thermal inflation is realized after the primordial inflation. When such a mini-inflation 
occurs, the time when the fluctuations of the reference scale in the present observations 
exit the horizon during inflation is shifted toward the end of inflation. This, in most cases, 
means that the spectral index is more tilted. Thus even if an inflation model is consistent 
with current observations in the standard case, some models can be excluded when the 
thermal inflation occurs. We have discussed that this situation can be relaxed by adopting 
the inhomogeneous end of thermal inflation since fluctuations from the mechanism are 
generally more scale- invariant as seen from Eq. f|T7|l . We have explicitly worked on some 
simple inflation models, such as chaotic inflation, new inflation and hybrid inflation. Then 
we find that some excluded models can be relaxed by using this mechanism proposed in 
this paper. 

As a final comment, we mention baryogenesis in a scenario with the thermal inflation. 
Since the thermal inflation dilutes the preexisting baryon asymmetry, one needs to gener- 
ate the baryon number after thermal inflation. Some authors have discussed Affleck-Dine 
baryogenesis after thermal inflation [31]. It may also be interesting to investigate observa- 
tional signature in such model with the inhomogeneous end of thermal inflation scenario, 
for example, baryon isocurvature fluctuations. This would be a subject of the future study. 

Acknowledgments 

This work is supported by Grant-in- Aid for Scientiflc research from the Ministry of Educa- 
tion, Science, Sports, and Culture, Japan, No. 14102004 (M.K.) and No. 19740145 (T.T.), 
and also by World Premier International Research Center Initiative, MEXT, Japan. S.Y. 
is supported in part by Grant-in-Aid for Scientiflc Research on Priority Areas No. 467 
"Probing the Dark Energy through an Extremely Wide and Deep Survey with Subaru 
Telescope" . He also acknowledges the support from the Grant-in-Aid for the Global COE 



12 



Program "Quest for Fundamental Principles in the Universe: from Particles to the Solar 
System and the Cosmos" from MEXT, Japan. 



References 

[1] K. Yamamoto, Phys. Lett. B 168, 341 (1986). 

[2] P. Binetruy and M. K. Gaillard, Phys. Rev. D 34, 3069 (1986). 

[3] G. Lazarides, C. Panagiotakopoulos and Q. Shafi, Phys. Rev. Lett. 56, 557 (1986); 
G. Lazarides, C. Panagiotakopoulos and Q. Shafi, Phys. Rev. Lett. 58, 1707 (1987); 
G. Lazarides and Q. Shafi, Nucl. Phys. B 392, 61 (1993). 

[4] D. H. Lyth and E . D. Stewart, Phys. Rev. Lett. 75, 201 (1995) 
jar Xiv:hep-ph/95024l7| . 

[5] D. H. Lyth and E. D. Stewart, Phys. Rev. D 53, 1784 (1996) |arXiv:hep-ph/95 10204] . 

[6] T. Barreiro, E. J. Copeland, D. H. Lyth and T. Prokopec, Phys. Rev. D 54, 1379 
(1996) [arXiv:hep-ph/9602263j . 

[7] T. Asaka and M. Kawasaki, Phys. Rev. D 60, 123509 (1999) |arXiv:hep-ph/9905467| . 

[8] G. D. Coughlan, W. Fischler, E. W. Kolb, S. Raby and G. G. Ross, Phys. Lett. B 
131, 59 (1983). 

[9] T. Banks, D. B. Kaplan and A. E. Nelson, Phys. Rev. D 49, 779 (1994) 
|arXiv:hep-ph/9308292 |. 

[10] B. de Carlos, J. A. Casas, F. Quevedo and E. Roulet, Phys. Lett. B 318, 447 (1993) 
|arXiv:hep-ph/9308325 |. 

[11] E. Komatsu et al. [WMAP Collaboration], Astrophys. J. Suppl. 180, 330 (2009) 
[arXiv:0803.0547l [astro-ph]]. 

[12] K. Enqvist and M. S. Sloth, Nucl. Phys. B 626, 395 (2002) [arXiv:hep- ph/0109214]. 



[13] D. H. Lyth and D. Wands, Phys. Lett. B 524, 5 (2002) [arXivihep-ph/Ol 1000 2 1. 

[14] T. Moroi a nd T. Ta kahashi, Phys. Lett. B 522, 215 (2001) [Erratum-ibid. B 539, 303 
(2002)] [arXiv:hep-ph/0 110096] . 

[15] G. Dvah, A. Gruzinov and M. Zaldarriaga, Phys. Rev. D 69, 023505 (2004) 
|arXiv:astro-ph/0303591j . 



[16] L. Kofman, arXiv:astro-ph/0303614 



13 



[17] E. W. Kolb, A. Riotto and A. Vallinotto, Phys. Rev. D 71, 043513 (2005) 
|arXiv:astro-ph/041Q546| . 

[18] D. H. Lyth, JCAP 0511, 006 (2005) farXiv: astro-ph/05 10443| . 

[19] L. Alabidi and D. Lyth, JCAP 0608, 006 (2006) |arXiv:astro-ph/0604569| . 

[20] M. Sasaki, Prog. Theor. Phys. 120, 159 (2008) jarXiv:0805.0974l [astro-ph]]. 

[21] K. M. Smith, L. Senatore and M. Zaldarriaea. larXiv:0901.2572l [astro-ph]. 

[22] N. Bartolo, S. Matarrese and A. Riotto, Phys. Rev. D 69, 043503 (2004) 
|arXiv:hep-ph7 0309033|; K. Enqvist and S. Nurmi, JCAP 0510, 013 (2005) 
arXiv:astro-ph /0508573 1 ; K. A. Mahk and D. H. Lyth, JCAP 0609, 008 (2006) 
arXiv:astro-ph/0604387|; M. Sasaki, J. Va hviita and D. Wands, Phys. Rev. D 
74, 103003 (2006) [arXiv:astro-ph/0607627]; Q. G. Huang, ■arXiv:08 01.0467 [hep- 
th]; K. Ichikawa, T. Suyama, T. Takahashi and M. Yamaguchi, Phys. Rev. D 78, 
023513 (2008) j arXiv:0802.4T38l [astro-ph]]; K. Enqvist and T. Takahashi, JCAP 
0809, 012 (2008) [ arXiv:0807.3069l [astro-ph]]: Q. G. Huang, JCAP 0811, 005 (2008) 
(arXiv:0808. 17931 [hep-th]]; Q. G. Huang and Y. Wang, JCAP 0809, 025 (2008) 
|arXiv:0808.TT68l [hep-th]]: T. Moroi and T. Takahashi, Phys. Lett. B 671, 339 (2009) 
|arXiv:08 10.0T89' [hep-ph]]; M. Kawasaki, K. Nakayama and F. Takahashi, JCAP 
0901 (2009) 026 |arXiv:0810.1585] [hep-ph]]; P. Chingangbam and Q. G. Huang, 
JCAP 0904, 031 (2009) |arXiv:0902.2"6l9l [astro-ph. CO]]: K. Enqvist and T. Taka- 
hashi, larXiv:0909.5362l [astro-ph. CO]. 

[23] M. Zaldarriaga, Phys. Rev. D 69, 043508 (2004) |ar Xiv:astro-ph/0 306006|; T. Suyama 
and M. Yamaguchi, Phys. Rev. D 77, 023505 (2008) |arXiv:0709.2"545l [astro-ph]]; 
K. Ichikawa, T. Suyama, T. Takahashi and M. Yamaguchi, Phys. Rev. D 78, 063545 
(2008) [arX iv:0807 .3988 [astro-ph]]; T. Takahashi, M. Yamaguchi and S. Yokoyama, 
farXiv: 0907.30521 [astro-ph. CO]. 

[24] T. Matsuda, Class. Quant. Grav. 26, 145011 (2009) |arXiv:0902.4283l [hep-ph]]. 

[25] A. A. Starobinsky, JETP Lett. 42 (1985) 152 [Pisma Zh. Eksp. Teor. Fiz. 42 
(1985) 124]; M . Sasaki and E. D. Stewart, Prog. Theor. Phys. 95, 71 (1996). 
jarXiv:astro-ph/9507001|; M. Sasaki and T. Tanaka, Prog. Theor. Phys. 99, 763 
(1998). [arXiv:gr-qc/9801017] . 

[26] D. H. Lyth, K. A Mahk and M. Sasaki, JCAP 0505, 004 (2005) 
|arXiv:astro-ph/0411220| . 

[27] D. H. Lyth and Y Rodriguez, Phys. Rev. Lett. 95, 121302 (2005) 
|arXiv:astro-ph/0504045]; L. Alabidi and D. H. Lyth, JCAP 0605, 016 (2006) 
arXiv:astro-ph/0510441|; C. T. Byrnes, M. Sasaki and D. Wands, Phys. Rev. D 
74, 123519 (2006) [arXiv:astro-ph/0611075] . 



14 



[28] D. Langlois and F. Vernizzi, Phys. Rev. D 70, 063522 (2004) 
|arXiv:astro-ph / 0403258 1 ; T. Moroi, T. Takahashi and Y. Toyoda, Phys. Rev. 
D 72, 023502 (2005) |arXiv:hep-ph/050100 7]; T. Mo roi and T. Takahashi, Phys. 
Rev. D 72, 023505 (2005) [arXiv:a stro-ph/0505339l ; K. Ichikawa, T. Suyama, 
T. Takahashi and M. Yamaguchi, in Ref. [22]; K. Ichikawa, T. Suyama, T. Takahashi 
and M. Yamaguchi, in Ref. [23]. 

[29] E. D. Stewart, Phys. Lett. B 345, 414 (1995) |arXiv:astro-ph/9407040l . 

[30] G. Lazarides and C. Panagiotakopoulos, Phys. Rev. D 52, R559 (1995) 
|arXiv:hi^ph/9506325j . 

[31] E. D. Stewart, M. Ka wasaki and T. Yanagida, Phys. Rev. D 54, 6032 (1996) 
|arXiv:h ep-ph/960332 4|; D. h. Jeon g, K. Kadota, W. I. Park and E. D. Stewart, JHEP 
0411, 046 (2004) |arXiv:hep-ph/0406136|; M. Kawasaki and K. Nakayama, Phys. Rev. 
D 74, 123508 (2006) |arXiv:hep-ph/0608335| ; G. N. Felder, H. Kim, W. I. Park and 
E. D. Stewart, JCAP 0706, 005 (2007) |arX iv:hep-ph/0703275| ; S. Kim, W. I. Park 
and E. D. Stewart, JHEP 0901, 015 (2009) |arXiv:0807.3607l [hep-ph]]; K. Choi, 
K. S. Jeong, W. I. Park and C. S. Shin. larXiv:0908.2154l [hep-ph]. 



15 



